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Abstract 

Neutrino non-standard interactions (NSI) were investigated earlier in the solar case 
and were shown to reduce the tensions between the data and the large mixing angle 
solution predictions. We extend the previous framework to the supernova and evaluate 
the appearance probabilities for neutrinos and antineutrinos as a function of their 
energy after leaving the collapsing star with and without NSI. For normal hierarchy the 
probability for electron neutrinos and antineutrinos at low energy (E < 0.8 — 0.9MeV) 
is substantially increased with respect to the non-NSI case and joins its value for inverse 
hierarchy which is constant with energy. Also for inverse hierarchy the NSI and non- 
NSI probabilities are the same for each neutrino and antineutrino species. Although 
detection in such a low energy range remains at present an experimental challenge, it 
will become a visible trace of NSI with normal hierarchy if they exist. On the other 
hand the neutrino decay probability into an antineutrino and a majoron, an effect 
previously shown to be induced by dense matter, is, as in the case of the sun, too 
small to be observed as a direct consequence of NSI. 
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1 Introduction 



Non-standard neutrino oscillations (NSI) [IHH] have since long been applied to solar neutri- 
nos [UEJIlQ] in an attempt to understand the origin of their apparent deficit. More recently 
interest in NSI has been revived [TT|[T2] with the purpose of solving the possible inconsisten- 
cies between the LMA solution to the solar neutrino problem [13j[TI] an d the data [T5T4T8] . 
The most remarkable of these inconsistencies is the absence of any experimental evidence 
for an upturn in the LMA survival probability in the intermediate energy range of solar 
neutrinos. It has been shown [TTJ that NSI can lead to a flat probability in this sector 
and hence a flat electron energy spectrum as observed in both SuperKamiokande [T3] and 
SNO [17] experiments, with the interesting consequence of a possible neutrino decay into 
antineutrinos and majorons in dense matter. However, the predicted antineutrino flux is 
too small to be observed and the only expected signature of the effect is the flatness of the 
electron spectrum. Extension to supernova neutrinos, which will be done in the present 
paper, will provide us further information on possible experimental NSI signatures. Our 
approach to NSI, previously developed in [TT] assumes extra contributions to the vertices 
v a Vp and v a e and differs therefore from the neutrino-neutrino interaction considered in [7J. 

The supernova dynamics has been extensively studied and for details we may refer the 
reader to the review [19]. Here we will just highlight its relevant aspects for our purposes. 
As is well known when a star of mass M > 8M has burned all its fuel, an onion like 
structure is formed with the lighter elements in the outer layers and the heavier ones inside. 
For massive enough stars M > 11M Q an innermost iron core is formed. Equilibrium of the 
core is disrupted by photodissociation of the heavier elements producing alpha particles and 
neutrons. The resulting free electrons are in turn captured by protons and nuclei, producing 
electron neutrinos which escape, 



When the Coulomb pressure of the electrons becomes insufficient to sustain the core against 
its own gravity, collapse is initiated until the core becomes no longer transparent to neutrinos 
and reaches nuclear density (p ~ 3 x 10 u g cm -3 ), a process which takes only about 10 
ms [20] . The neutronization process ([I]) ceases, a state of hydrostatic equilibrium is reached 
with the core forming a proto-neutron star of radius around 10 km and a temperature of 
35-40 MeV. From this stage onwards only thermal neutrinos are emitted. They are produced 
from nucleon-nucleon bremsstrahlung 
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and electron nucleon bremsstrahlung 

e~ N e~ N v v. (5) 

So according to core collapse supernova dynamics, in the thermal phase neutrino emission 
proceeds through all flavour channels and is accompanied by antineutrinos, whereas in the 
initial neutronization phase only electron neutrinos are expected. 

The events observed in SN1987a were all interpreted as antineutrino ones [21] due not 
only to the dominance of the antineutrino cross section relative to the neutrino one but also 
possibly to the unavailability of a lower energy threshold, as shall be seen in the present 
paper. Whether or not any neutrinos reached the Earth is an open question. In ref. [11] 
it was shown that neutrino decay in dense matter into antineutrinos and a majoron is a 
necessary consequence of NSI. However its effect through the detection of antineutrinos 
in the neutronization phase where only z/ e 's are produced is not possible, owing to the 
smallness of the antineutrino production probability through decay even in the supernova. 
In fact the high density range of the neutrino trajectory is not long enough to produce a 
visible antineutrino flux originated from this source. 

In this paper we extend for supernova neutrinos the NSI framework that was introduced 
in the solar neutrino case [llj where it was shown to improve the fits to experiment. In 
section 2 we review its general features. We recall that a flat SuperKamiokande and SNO 
electron spectrum necessarily require imaginary diagonal entries in the NSI Hamiltonian. 
Their real parts, along with the off diagonal entries whether real or imaginary, may be 
arbitrary in the sense that they do not induce any change in the standard LMA probability. 
As a consequence neutrino decay in dense matter into an antineutrino and a majoron arises. 
Section 3 is the main part of our work. Here we extend NSI to the interactions of the 
supernova neutrinos. We evaluate the survival and conversion probabilities with and without 
NSI as well as the probability for neutrino decay. For normal mass hierarchy in the absence 
of NSI the electron neutrino survival probability turns out to be small, most of the i/ e 's 
having been converted through standard oscillations to z/ M , i/ T 's. Furthermore it is found 
that as in the case of the sun, the decay probability is extremely small for antineutrinos 
from neutrino decay to ever be observed. Since in the neutronization phase only u e 's are 
produced, no charged current signal is expected to be seen at this initial stage in the absence 
of NSI. On the other hand in the NSI case a sizeable electron neutrino flux may appear at 
low energies in the neutronization phase which may be detected through the charged current 
with improved low energy detectors. At present this remains a challenge but may be reached 
in the not too distant future. For inverse hierarchy all neutrino fluxes are comparable in 
the whole energy range. On the other hand in the thermalization phase neutrinos (i/ a ) and 
antineutrinos (p a ) of all kinds are produced so that the initial state is assumed to consist 
of u a , Va s in equal proportions. Owing to the large number of oscillations, the information 
from the initial state is essentially lost, hence the final probability distributions are the 
same as in neutronization where only v<?s are present initially. The only difference is in 
this case the obvious presence of antineutrinos with the same probability distribution as 
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neutrinos due to CPT invariance. Finally in section 4 we comment on our results and draw 
our conclusions. 



2 The NSI Hamiltonian 



In this section we highlight and discuss the main steps of the analysis leading to the Hamil- 
tonian for propagation in dense matter (see also [H]). Its results will be used in section 3 for 
the evaluation of the probabilities for neutrino survival and conversion and the probability 
for antineutrino production. 

Assuming only standard interactions (SI), the potential for electron neutrinos traveling 
through the sun and supernova is given by 



V = G F V2N e 1 1 - 
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where N e , N n are the electron and neutron density and V c = Gpy/2N e , V n = —Gp/y/2N n . 
For NSI we assume that the u a interaction potential on electrons (a = e, /x, r) involves both 
the charged and neutral currents (CC and NC), while on quarks it involves only NC. NSI 
give rise to possible lepton flavour violation. Denoting by e^l the NSI factor that multiplies 
each diagram associated to neutrino propagation in matter we have 
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In the following we will assume (e^)cc — ( e ap)NC- Hence flavour change may occur without 
a vacuum mixing angle [13JH1] or a magnetic moment [22], being induced only by the off 
diagonal entries of this matrix (a ^ f3). So with SI and NSI the matter Hamiltonian in the 
flavour basis is the sum of eqs.([6]) and (J7J) 
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As is well known, owing to the large neutrino density, collective effects can play an important 
role in the supernova [23], HU, [22] @, providing an additional contribution to the matrix (JHJ). 
They occur up to a few 100 km, whereas MSW oscillations occur typically at larger distances 
so that MSW effects factorize and can be included separately [27]. For this reason the NSI 
effects, intrinsically associated to MSW, must also be included separately. We will perform 



*For a review on collective oscillations see 
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our calculation starting from a region around 400-500 km where collective oscillations have 
already taken place. Hence their net effect to our approach amounts to the well known 
spectral swap-split of the neutrino and antineutrino energy spectra [24], [25], [28]. However 
it was recently shown that matter completely suppresses collective oscillations up to 200 
ms [29J after bounce. In view of these results it appears that collective effects can be 
ignored at early times in a supernova. Our results, derived in the next section, are therefore 
expected to be valid for the whole neutronization phase and part of the subsequent thermal 
phase. 

The investigation performed in ref. [11] shows that no off diagonal entry in matrix 
(jSJ) whether real or imaginary, can change the LM A probability, hence the rates and the 
corresponding SuperKamiokande and SNO electron spectra. In fact only imaginary diagonal 
couplings lead to a change in Plma- As discussed in ref. [IT] , this implies the instability of 
neutrinos in matter and their decay into antineutrinos of all species along with majorons. 
Requiring the convenient change in Plma, namely the one that leads to a flat electron 
spectrum in SuperKamiokande and SNO, thus allowing for imaginary diagonal couplings, it 
was shown that the simplest choice of parameters is 
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with vanishing off diagonal entries and £ = 3.5 x 10 4 . Here 

1 - 8 2 „ JV„ 1 2 2 „ JV„ 

x e = -- + 2sm 6 W , x u = 1 - -sin W + — , x d = -- + -sm 6 W - — . (10) 

The three matrices in the right hand side of eq. ([9]) relate to the neutrino interaction with 
electrons, u-quarks and d-quarks respectively. Each diagonal entry refers to the u e , u^, v T 
contribution to its own interaction, hence its decay. So for instance, for u e the NSI (decay) 
potential is 

(v ee ) NS i = l -eG F y/2N e (x e + 1 + x u - x d ) . (11) 
In this way we obtain the three interaction (decay) potentials 

Mnsi = ^^2(3.5 x 10- 4 )iV e f 1 - -sinH w + j^j (12) 

(v^) NS i = iG F V2{3.5 x 10- 4 )iV e ^-1 + ~sin 2 6 w - (13) 
Mnsi = ^G F V2(3.5 x 10' 4 )iV e ^1 - Uin 2 6 w + . (14) 
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The full Hamiltonian including the vacuum part and referred to the mass basis is now 
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with Hm given by (jSJ) with the replacements ( IT2l) . ( {TBI . ( Fl4l) and arbitrary off diagonal 
entries. In eq. ( fl5]) Am^, Am| x are the solar and atmospheric mass differences, U is the 
PMNS matrix [30] defined with the standard parameterisation [31] and we use the central 
value for sin8\^ = 0.13 from ref. [13]. The negative imaginary parts of the eigenvalues of 
( IT5l) are the mass eigenstate decay rates Tj which will be used in the next section for the 
evaluation of the probabilities. As for the supernova parameters, numerical simulations [32] 
yield the electron number density and supernova density profiles which in our period of 
interest (the initial 200 ms) are well approximated by 
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for r > 10 km with po ~ 10 g cm and which will be used in the next section. 



3 Rates, couplings and probabilities 

3.1 Rates and neutrino majoron couplings 

The v e flux from neutronization is in fact a linear combination of the three mass eigenstates 
Vi displayed in figfjja) for neutrino energy Eq = HMeV and normal hierarchy. Our first 
purpose in this section is to evaluate the decay rate of the NSI process [33H38] 

»i ^&j + X (18) 
where x denotes the majoron. This rate satisfies [TT] . [37] 

§ r t^^Mr)-^r)\, SI F(r,E 0) (19) 
i=i 

where are the neutrino majoron couplings [37], [SS] and the interaction potentials satisfy 
in the mass basis Vi = — Uj. Here E is the initial neutrino energy, Ef is the antineutrino 
energy which in a first approximation we assume to take values in the interval (0, E ), since 
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the energy Eq is shared by the final neutrino and the majoron. The quantity F(Eq, r) is the 
Fermi factor [37J 

F ( E ^=( 1 - e ( Eo -')/T + 1 ) > ( 2 °) 

which reflects the fact that inside the supernova some of the states have already been 
occupied by neutrinos. In the inner core (Rinner — 10 km) the chemical potential for v e (fi Ue ) 
is around 200 MeV and the temperature T = 35 MeV. In the outer core (Ri nner — 15 km) 
the temperature drops abruptly T = 2 MeV, the density falls from 5 x 10 14 g cm" 3 to 
5 x 10 13 g cm -3 and we may set \i Vts = 0. Hence in the following we will omit the factor 
F(E ,r). 

The three decay rates Tj are the imaginary parts of the NSI Hamiltonian eigenvalues. 
They are represented in figJT]^) for Eq = HMeV. As in the solar case, only 1^ is negative, 
so only the state u 3 is unstable, allowing for the decay into either Pi, z/ 2 or P 3 , thus generating 
antineutrinos of the three flavours. One could obtain an alternative expression for Tj(r, Eq) 
through the integration of equation (j!9p over Ef 

r ' = Cw f dE '=tM Mr) ~ §i{r)l - (21) 

The vanishing lower limit used in this integration is as referred to above only an approxima- 
tion, since the majoron obtains a tiny effective mass in matter vr? e ^ ~ \g\ 2 N u /q [37J. In this 
way a slight dependence of the rate Tj on E arises which is consistent with the numerical 
evaluation of the imaginary part of the NSI Hamiltonian eigenvalues. 

The parameters e a p were fixed earlier (see section 2) by the fittings to the solar neutrino 
data [H] and so the rates Ti{r, Eq), numerically evaluated as the imaginary parts of the 
Hamiltonian eigenvalues, are also fixed. Moreover the values of gij to which the rates 
correspond are so far unknown and only upper bounds exist in the literature [38J. In the 
following we will use the strictest one quoted, namely 

^|£ eQ | 2 <5.5xlO- 6 , (22) 

a 

conveniently expressed in the mass basis. Quantities Tj will now be used in the evaluation 
of the probabilities. 



3.2 Probability densities for neutrino survival, decay and antineu- 
trino appearance 

We denote by P v Sjt -^o) the ^ survival probability in the mass basis for neutrino energy Eq 
at a distance r from the star centre, obtained from integration of the Schrodinger equation 
with the Hamiltonian (j!5p . and by <p Ue (E ) the initial normalized neutrino spectral flux [20] 

1 9$ 
^ {Eo) ~ HEo)dE-Q ■ 
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Given these definitions the quantity 

<f) Ui (r,E ) = P Ui (r,E )<p Ue (E ) (23) 

is the normalized spectral flux of V{ mass eigenstates with energy Eq that remain in the 
beam after traveling a distance r. Hence 

M = f M r, - ^) »^ ^ (24) 

is the probability per unit star radius for this mass eigenstate to have disappeared from the 
flux. 

Replacing now ( Fl9|) in ff24|) one obtains the probability for the mass eigenstate disap- 
pearance 

Pu T {E Q ) = f^£- f RS \ Vl (r) - vMxsiMr, E )(l - e^>) ^ldE f dr (25) 

Quantities Ri and Rs denote the neutrino sphere and the star radii respectively. We note 
that in eq.f )25|) we have the simple sum of probabilities, since each eigenstate z/ i; as it decays, 
can give rise to just one antineutrino flavour which is a linear combination of the three mass 
eigenstates Uj. Performing the integration in Ef, (T2"5"j) can be simplified to 

P ur (E )= fj^f RS \ Vt ( r )-v 3 {r)\ NSI <j> Vi {r,E Q ){l - e r ^ r ) dr . (26) 

jr[ 167r Jr, 

In other words, given a flux of s with an energy in the interval [Eq, E + dEo], the quantity 
P u m[Eo)dEo is the fraction of these neutrinos which has decayed into antineutrinos with an 
energy in the interval (0, Eq + dEo) after traversing the star. The fraction of Vi s that 
remains in the beam after leaving the star is P Ui (R s , E )dE with P u .(r,E ) as defined in 
the beginning of this section (see (123]) ). 

As a reminder we note that the normalization of P Vi follows from 

P Vi = \< Vi\v e > | 2 = \viuljVjl 2 = uljUekVjVkViVi = u* ej u ek 5ijS ik = u* ei u ei (27) 

(no sum over i) with the orthogonality of JJ PMNS ensuring ^ P Vi = 1. Moreover the 
integration of the equation of motion proceeds in the mass basis, so that the normalization 
condition is forced on P u . . The normalization of P Va also follows in a similar way 

Pv a = | < V a \v e > | 2 = \u* a i u ejViVj\ 2 = U* ai U* ej U ak U em PiVjV k V m 

(no sum over a) with the orthogonality of JJ PMNS implying ^ a P„ a = 1. In (1271) and (1281) 
we neglected the plane wave propagation phases. Finally the flavour probability is obtained 
from 

P Va = \u ai < Vi\v e > | 2 (29) 
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and is shown in figl2]for standard (e = 0) and non-standard interactions (e = 3.5 x 10 4 ). 

For the sake of the following discussion we recall that from supernova theory only elec- 
tron neutrinos are produced in the initial neutronization phase. From figj2] with normal 
hierarchy (panel (a)) it is seen that, in the absence of NSI, electron neutrinos can hardly 
be detected, as the survival probability is around 0.02 and practically constant with energy. 
Further, v^s and v^s produced from oscillations, although having a much larger production 
probability (~ 0.5), can only be detected through neutral current reactions, so there will 
be no way to distinguish them at present from other neutrinos or antineutrinos. The major 
difference in the NSI case for normal hierarchy as seen from figj2^, is the appearance of a 
visible flux of electron neutrinos in the low energy region (Eq < 0.8 — 0.9MeV). This can 
be pinpointed in the neutronization phase through the charged current, provided enough 
experimental capability is developed in the future to reduce the low energy threshold. On 
the other hand for inverse hierarchy (figfjb) no difference appears between NSI and non- 
NSI: the sizable value of the v e probability (~ 0.31) is the same at low energy as in the 
normal hierarchy case and remains constant with energy, so v e detection is experimentally 
accessible. However the NSI and non-NSI probability values are the same for each neutrino 
species and also remain in each case constant with energy. So for inverse hierarchy NSI and 
non-NSI appear indistinguishable. 

We next investigate the possibility for electron antineutrinos to be detected in the ini- 
tial neutronization phase. As previously referred, only electron neutrinos are produced at 
this stage in a 10 ms pulse. Their decay in matter into electron antineutrinos in a signif- 
icant number would provide a clear signature of NSI in this initial pulse. Antineutrinos 
with energy Ef are produced from electron neutrino decay with energy E in the interval 
(Ef, E 0max }. Their appearance probability density is in the mass basis 

, .3 \n--\ 2 f Rs r E o m ax , e\_/t\ 

^(£/)=U nr-/ K^-vMnsi <p u Ar,E )(l-e r ^ E °>)^-^dE dr. (30) 

87T J Ri J Ef £j 

Contrary to eq. (T2"6l) . where each neutrino could decay into one single antineutrino, each 
antineutrino can be now produced from more than one neutrino, hence the reason to consider 
the union of events in eq. dsn])!. Assuming CPT invariance the flavour probability for v e is 
given by 

Pp e = \u ei \ 2 P Pt • (31) 

Using the bound (|22|) as a common value for all neutrino majoron couplings involved, we 
obtain the result displayed in figj3]for P Pe which is obviously too small for v^s to be observed. 
If one considers the v e decay into v T one may relax ( 1221) and use instead [38] 

^|(?™| 2 <5.5xl0- 2 (32) 



§The probability for the union of three independent events (Ai,A 2 , ^3) is given by the well known rule 
P(A) = P(Ai) + P(A 2 ) + P(A 3 ) - P(A 1 )P(A 2 ) - P(A 1 )P(A 3 ) - P(A 2 )P(A 3 ) + P(A 1 )P(A 2 )P(A 3 ). 



S 



as a common value for all couplings. In this case the uppermost value of the probability is 
raised by a factor O(10 4 ) which is still too small for the effect to be observed and moreover 
i>eS, if produced from v T decay, would appear as a higher order effect. On the other hand, as 
pointed above, the v e probability from oscillations in the absence of NSI is seen to be rather 
small in normal hierarchy (figJ2K). The prospects for its detection and hence a charged 
current signal in the neutronization phase will crucially depend on the detector size and 
supernova distance. With NSI, it may become clearer through u e d — > ppe~ or an increased 
ve~ ve~ scattering event rate, however at low energy (E < 0.8 — 0.9MeV), which is 
still an experimental challenge. For inverse hierarchy, as also pointed above, the situation 
is much different: the v e signal appears louder and clearer (see figj2b). 

In the subsequent thermalization phase v e 's and vx^x's with X = fj,,r are produced 
initially along with v^s. All kinds of neutrinos and antineutrinos will arise from these 
through NSI. The appearance probability for antineutrinos from neutrino decay may increase 
substantially, since the bound ( 1321) must now obviously be taken into account instead of (1221) . 
Again this is not expected to be enough for the effect to be observable (see figEJ) even for z/ e 
appearance, despite the higher rate involved in its detection. As for the other probabilities 
from oscillations and NSI, they must be evaluated from the union of events as in (1301 . since 
each final neutrino or antineutrino is produced simultaneously from a number of different 
initial ones. Evaluating these probabilities taking into account their normalization, it turns 
out that the contribution from the extra initial neutrinos and antineutrinos does not change 
their value nor energy distribution relative to the neutronization case when initially only i/ e 's 
are present. In fact the information from the initial state is lost due to the exceedingly large 
number of oscillations that the neutrinos undergo during propagation: only the propagation 
physics which depends on the interaction potentials is relevant here. Thus the same figf5] 
applies for thermalization both for normal and inverse hierarchies. The only notorious 
characteristic to tell NSI from non-NSI is, as in neutronization and normal hierarchy, the 
comparatively large probability for u e , v e at low energy (Eq < 0.8 — 0.9MeV) with the same 
energy profile. Such low energy raises however an experimental challenge for detection. For 
inverse hierarchy u e 's and v<?s remain equally abundant both for Eq < 0.8 — 0.9MeV^ and 
larger, therefore their detection is acessible, although NSI and non-NSI are indistinguishable 
in this case (see figf2b). In particular v^s provide a clear signal through the reaction 
v e p — > ne + whose cross section is O(10 2 ) larger than for scattering with electrons. Since 
thermalization is a much longer process than neutronization (>10 s), a larger accumulation 
of events is possible in this phase. 

For the detection and measurement of the v T , v T individual energy spectra which 

can only be traced via neutral currents, an interesting proposal was presented some time 
ago [39] and recently revived [40J. It is based on the v p — > v p scattering reaction which 
can be observed in scintillator detectors (e.g. Borexino, SNO+, KamLAND) through their 
adequate preparation. This is a neutral current process with a cross section about O(10 2 ) 
larger than neutrino electron scattering at supernova neutrino energies. For the NSI scenario 
expound in the present paper this technique will be particularly useful, since it appears to 
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be possible to clearly distinguish between normal and inverted hierarchies. In fact it suffices 
to note that for normal hierarchy the above mentioned neutrinos arrive copiously on Earth 
in comparison with the more rare u e 's and u e 's whereas for inverse hierarchy all species arrive 
in comparable numbers (see figEk, b). 

4 Summary and conclusions 

We have extended to the supernova the previously developed model for neutrino NSI in 
the sun introduced earlier to remove the tension between the LMA predictions and the 
experimental signatures of solar neutrinos, especially the absence of an upturn in the Su- 
perKamiokande event rate. Improving the data fittings in the solar case implies neutrino 
decay in dense matter into antineutrino and a majoron, hence the motivation to investigate 
the consequences of the model for the supernova. 

In the present paper we found however that, although the matter density in supernova 
is much larger than in the sun, the extension of neutrino trajectory in the very high density 
medium is too short to imply a significant neutrino decay into antineutrino and a majoron 
and the corresponding appearance probability is insignificant. In the initial and short neu- 
tronization phase (10 ms) where only neutrinos are produced through the reaction (CQ) no 
antineutrinos are expected experimentally whether or not NSI applies. The important NSI 
trace is the v e appearance probability which increases from 0.02 to 0.31 at low energies 
(Eq < 0.8 — 0.9MeV) for normal hierarchy, while for inverse hierarchy NSI and non-NSI 
cannot be distinguished. In this case the v e probability remains at 0.31 regardless of the 
energy. In the neutronization (deleptonization) phase u e 's are the only states that can in- 
duce charged current interactions, so they can be singled out either through an increased 
ve~ — > ve~ scattering event rate or v e d — > ppe~ . The remainder (v^s, v^s) inducing only 
neutral currents, cannot be distinguished from each other nor from i/ e 's. Detecting these 
UeS remains however an experimental challenge at present in normal hierarchy, but not so 
for inverse hierarchy, as they appear more copiously at higher energies. In the subsequent 
and longer thermalization phase, the extra neutrino and antineutrino states (u e and ux, 
with X = /i, r) that are produced through processes (2)-(5) cannot change the appearance 
probabilities relative to the neutronization phase. Hence detecting these z/ e 's and i> e 's is, 
again, an experimental challenge in normal hierarchy: as for the neutronization phase v e 's 
can be detected through a major event rate increase originated from the charged current 
in ve~ — > ve~ scattering or through the reaction u e d — > ppe~ , while v^s through the clear 
signal v e p — > ne + . The remaining neutrinos and antineutrinos can only be detected through 
the neutral current and so cannot be distinguished from u e 's and u e 's, unless the interesting 
technique proposed in [39], [40] is developed. If and when this advancement succeeds, it may 
be possible within the present scenario to tell normal from inverse hierarchy. As regards 
collective oscillations, their effect in our analysis amounts to the modification of the neutrino 
and antineutrino spectral fluxes. Collective effects are however probably suppressed up to 
0.2 seconds after bounce. The results obtained are therefore applicable to the neutronization 
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(deleptonization) phase and part of the subsequent thermal phase. 

To summarize, in the presence of NSI we expect a sizable flux of v^s and z/ e 's at all 
energies for inverse hierarchy and at low energy (E < 0.8 — 0.9MeV) for normal hierarchy. 
These fluxes are the same as for non-NSI. The clear distinction between NSI and non-NSI 
is possible only for normal hierarchy at low energy with a more intense flux of v e and 
v ei whose detection is at present an experimental challenge. Hence in the absence of NSI 
the chances for observation in normal hierarchy of a charged current signal do not appear 
much favourable at present, but they will of course mainly depend on the detector size and 
supernova distance. Other neutrinos and antineutrinos are in contrast abundantly present, 
however they can only be detected through the neutral current. As in the case of the sun 
the antineutrino appearance probability from NSI neutrino decay is in all cases too small 
for antineutrinos to be detected from this origin. 
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Figure 1: The real (a) and imaginary parts (b) of the neutrino mass matter eigenvalues for 
neutrino energy E Q = HMeV. The two resonances ('atmospheric' and 'solar') are clearly 
visible in panel (a) and the two positive imaginary parts are superimposed in panel (b). 

14 



0.6 



nr E non-NSI 

0.5 = 



E non-NSI 



i > 



0.4 p NS | 
E NSI 



0.3 = 



0.0 



i — i — i — i i i 1 1 1 1 — i — i — i i i 1 1 1 1 — n 



= NSI 
0.2 i- 

0.1 E- 

E non-NSI 



NH 



P v or P - 



(a) 



P„ orP 



I I H 



10 w 



10 c 



10' 



E v (eV) 



0.6 



0.5 =- 



- NSI and non-NSI 



0.4 — 



i > 
Cl 



z NSI and non-NSI 



0.3 — 



0.2 — 



0.1 =- 



0.0 



T 1 — I — I I I I I I 1 1 — I — I I I I I I 1 1 — q 



- NSI and non-NSI 



IH 



P„ orP 



P v or P - 



(b) 



J I I I I I I I I I zl 



1 tr- 



io 



10' 



Figure 2: The appearance probabilities P Ue , Pp e , P u ^, Pp M , P Ut , P Dt with and without NSI 
((a) normal hierarchy, (b) inverse hierarchy). For normal hierarchy and energy below 0.9 
MeV the NSI probability merges with its non-NSI value. 
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Figure 3: The antineutrino appearance probabilities from neutrino decay in matter through 
NSI. The upper bound |HP (the strictest) was used as a common value for the neutrino ma- 
joron couplings. Antineutrinos from neutrino decay in supernova matter cannot be detected 
even if the more conservative bound (CHj) is used. 
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